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Abstract 

> , 

' In this article, we study the perturbational method to construct the non-radially symmet- 

■ 

ric solutions of the compressible 2-component Camassa-Holm equations. In detail, we first 
fS) ' combine the substitutional method and the separation method to construct a new class of 

, analytical solutions for that system. In fact, we perturb the linear velocity: 

, , , u = c{t)x + b{t), (1) 

X 

u ■ 

5t 1 S'lid substitute it into the system. Then, by comparing the coefficients of the polynomial, we 

can deduce the functional differential equations involving {c{t),h{t), p^{Q,t)). Additionally, we 
could apply the Hubble's transformation 

- HI' 

to simplify the ordinary differential system involving (a(3f), p^(0, t)). After proving the 
global or local existences of the corresponding dynamical system, a new class of analytical 
solutions is shown. And the corresponding solutions in radial symmetry are also given. To 
determine that the solutions exist globally or blow up, we just use the qualitative properties 
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about the well-known Emden equation: 



dt' 



a(3t) 



T 



< 



a3 (34) 



(3) 



a(0) = ao > 0, d(0) = ai 



Our solutions obtained by the perturbational method, fully cover the previous known results 
in "M.W. Yuen, S elf- Svmilar Blowup Solutions to the 2-Component Camassa-Holm Equations, 
J. Math. Phys., 51 (2010) 093524, 14pp." by the separation method. 
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1 Introduction 

The 2-coniponent Camassa-Holm equations of shallow water system can be expressed by 



Here u = u{x, t) € R and p — p{x, t) > are the velocity and the density of fluid respectively. The 
constant cr is equal to 1 or — 1 . If cr = — 1 , the gravity acceleration points upwards [2] , [3] , [8] , [lO] 
and [9]. For a — 1, the researches regarding the corresponding models could be referred by [4], [6], 
[TU] and [S]. When p = 0, the system returns to the Camassa-Holm equation [T]. The searching 
of Camassa-Holm equation can capture breaking waves. Peaked traveling waves is a long-standing 
open problem |17) . 

In 2010, Yuen used the separation method to obtain a class of blowup or global solutions of 
the Camassa-Holm equations ^22^ and Degasperis-Procesi equations ^23^. In particular, for the 




(4) 



with 



m = u — o?u. 



'XX ■ 



(5) 
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integrable system of the Camassa-Holm equations with cr = 1, we have the global solutions: 

p{x,t) =max{^^^||2}^,0}, u{x,t) = 

- = 0' = «o > 0, a(0) = ai (6) 

where ry = ^^.^^^yg with s = 3t; ^ > and a > are arbitrary constants. |22j 

Meanwhile, the isentropic compressible Euler equations can be written in the following form: 



Pt+V • pu =0 
(pm)j+V • {pu u) + VP = 0. 



(7) 



As usual, p = p{x,t) and u = u{x,t) £ are the density and the velocity respectively with 
X = {xi, X2, xjv) G i?^. For some fixed K > 0, we have a 7-law on the pressure P — P{p), i.e. 



P{p) = Kp^ 



(8) 



with a constant 7 > 1. For solutions in radially symmetry: 



X X 

p{x,t) = p{r,t) and u{x,t) = —V{r,t) =: —V 

r r 



(9) 



where the radial r — J2i=i 



(10) 



the compressible Euler equations ([7]) become, 

pt + Vpr + pVr + = 

p{Vt+VVr)+VP = Q. 

Recently, there are some researches concerning the construction of solutions of the compressible 
Euler and Navier-Stokes equations by the substitutional method [12], [13]) [IB] and [Tl]. They 
assume that the velocity is linear: 

u{x,t)^c{t)x (11) 

and substitute it into the system to derive the dynamic system about the function c{t). Then 
they use the standard argument of phase diagram to drive the blowup or global existence of the 
ordinary differential equation involving c{t). 

On the other hand, the separation method can be governed to seek the radial symmetric solutions 
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by the functional form: 



, f^W)) , a{t) 



(12) 



(H, m. 0, ffu, m. m. m and nj) 

It is natural to consider the more general linear velocity: 

u{x,t) = c{t)x + h{t) 



(13) 



to construct new solutions. In this article, we can first combine the two conventional approaches 
(substitutional method and separation method) to derive the corresponding solutions for the sys- 
tem. In fact, the main theme of this article is to substitute the linear velocity (IT3)) into the 
Camassa-Holm equations ^ and compare the coefficient of the different polynomial degrees for 
deducing the functional differential equations involving (c(t), b{t),p'^{0, t)). Then, we can apply the 
Hubble 's transformation 

d(3t) 



c{t) 



a(3t) 



(14) 



with d(3i) :~ '^'^lif^ i to simplify the functional differential system involving {a{3t) , b{t) , {0 , t)) . 
After proving the local existences of the corresponding dynamical system, we can show the results 
below: 



Theorem 1 For the 2-component Camassa-Holm equations there exists a family of solutions, 



p\x,t) = max^ p2(o^^) _ l^r)+ib{tr^]^^z^y, 



era 3 (t) 



^a(3i) = a(0) = > 0, d(0) = ai 



a^(3t) 



&m + '-§^ib{t) + -|2L_&(t) = 0,6(0) = bo, 6(0) = b, 
i [pHO,t)] = pHO,t) ~ f [6(0 + 3b{t)^]x - p^(0,0) = 

where Oq, ai, bi, 62 and a are arbitrary constants. 



(15) 



We remark that the above solutions ((T5|) fully cover the previous known results (35] by the 
separation method by choosing 69 = 61 = 0. 
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2 Perturbational Method 

The proof for Theorem [1] requires the standard manipulation of algebraic computation only: 
Proof. The momentum equation becomes 

First, we perturb the velocity with this following functional form: 

u{x,t) = c{t)x + b{t) 

where b{t) and c{t) are the time functions determined later. 
As the velocity u (IT7)) is linear: 

it can be simplified to be 

Ut + 'iuUx + (Jppx = 

t{t)x + h{t) + i[c{t)x + b{t)]c{t) + -T-P^ = 

z ox 

f = -m + mxt)] m + 3c'{t)]x. 

Then, we take integration from [0,2;] to have: 
d 



a 



?ds = ~[b{t) + 3&(t)c(t)] I ds~ [c{t) + 3c^(t 



I [p^x^t) - p^m] = -m + 3b{t)c{t)]x - Mi^£!M^2 

p^{x,t) = p\Q,t) - -ib{t)+ib{t)c{t)]x - M±MMa;2^ 
a a 

On the other hand, for the 1-dimensional mass equation (j4])i, we obtain 

pt + [c{t)x + b{t)] px + pc(t) = 0. 
Here, we multiple p on both sides to have 
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After that, we can substitute equation ([M| into equation (pS)): 



2 vat 



9t 0- 



+ [c{t)x + b{t)] { --[b{t) + 3b{t)c{t)] - -[c{t) + 3c''{t)]x 
a (J 



c{t) 



p2(o, t) - -[h{t) + ih{t)c{t)]x - M±^£!Ma;2 



(27) 
(28) 
(29) 



^ I [P' (0, ^ )] + c(t)p^ (0, t ) - ^ [6(t ) + 36(t)c(t )] 



-if [fe(i)+36(t)c(t)] 



[6(t) + 3b{t)c{t)] 



b{t) 



'c{t) + 3c2(t)] - + 3&(t)c(t) 



(30) 



(31) 



-_L|[c(t) + 3c2(t)]-i[c(t) + 3c2(t)]c(t) 



c(t)[c(t)+3c^(t)] 



(32) 



By comparing the coefficients of the polynomial, we require the functional differential equations 
involving {c{t),b{t), p'^(Q,t)): 



A r^2 

dt 



p2 (0, t)] + 2c{t)p^ (0, t) ~ Hit) m + 36(t)c(t)] = 



I [6(0 + 36(t)c(t)] + 3c{t)[b{t) + 3b{t)c{t)] + b{t)[c{t) + 3c^{t)] = 
fjc{t) + 3c2(t)] + 4[c(t) + 3c2(t)]c(t) = 



(33) 



For details (existence, uniqueness and continuous dependence) about general functional differential 
equations, the interested reader may refer to the classical literatures [TT] and [TB]. 

For solving the above ordinary differential system (1331) . we initially solve equation ^31) -^ about 
the function c{t). Here we let the function c{t) be expressed with the Hubble's transformation: 



c{t) 



a{3t) 



(34) 



where d(3t) :— and the function a(3t) could be determined later. 



It is transformed to be 



d 3d{3t) 3a2(3i) 3a'^{3t) 3a(3t) 3d^{3t) 3d^{3t) d{3t) _ 

dt^^(3i)" ~ a^{3t) ^ a^{3t) ^ ^ ^^(3i)" " a'^{3t) ^ a^{3t)^^^~ ^ ' 
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d_ ( a(3t) \ 4a(3t) d(3t) _ „ 
dt \a(U) J a(3t) 0(3*) ^ 

a(0) = oq > 0, a(0) = ai, a(0) = a2 
3'a(3t) 3a(3i)a(3i) 4a(3t)a(3t) 

a{3t) a(3i)a(3i) 



= 



a{3t) 3a2(3i) 



0. 



Then, we multiple a^(3t) on both sides to have: 



a{3t) a (3t) + ' = 0. 



(36) 

(37) 
(38) 



(39) 



It can be reduced to the second-order Emden equation: 



a(3t) 



a3 (3t) 



a(0) = ao > 0, a(0) = ai 



(40) 



where ^ := aga2 is an arbitrary constant by choosing 02- 

We remark that the well-known Emden equation is well studied in astrophysics and mathematics. 

Next, for the second equation (IM])? about b{t) of the functional differential system, we could 
further simply it in terms of the known function a{3t): 



d r! , ^ „wxa(30i „a(3^)r,■/^ „wxa(30i 3a(3t) , , , 
OT a(3t) a(3t) a[3t) a[3t) 



a[3t) 



a(3i) _a2(3t) , 9a2(3i) , 3a(3t) 
9 — ■ , — 9 ' ' 



a{3t) "a^{3t) a^i3t) a{3t) 

12? 



b{t) = 



.d(3t) , 
'^(3ty' 

6(0) = 60, 6(0) = 61 



with the Emden equation (HUl) . 



We denote /i(t) = 6^ and /2(t) = 4^ to have 

^ ^ "^(S*) ^ ^ a3(3t) 

6(i) + /i(i)6(i)+/2(i)6(i) =0 
6(0) = 60, 6(0) = 61. 

Therefore, when the functions fi{t) and f2{t) are bounded, that is 



(41) 
(42) 

(43) 



(44) 



|/i(t)| <Fi and \f2{t)\ <F2 



(45) 
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with some constants -Fi and F2, provided that the functions ^^^^^ and a(3f) exist, the functions 
h{t) and h{t) can be guaranteed for existing by the comparison theorem of ordinary differential 
equations [T7] . 

Lastly, for the first equation (|33|) i . we denote H{t) — ^^^^ and G{t) — + 2ih{t) "jg^j ] in 

terms of functions ^|jy, a(3i), h{t) and h{t) provided that they exists, to solve 

^[p2(o,t)]+p2(o,t)i/(i)^G(i) 

(46) 

p2(0,0) = a2. 

The formula of the first-order ordinary differential equation (j46|) is 

9, ^ f!l ii(s)G(s)ds + k , ^ 

where 

Therefore, we have the density function from equation (j24p : 

= p2(0,i) - -W) + 36W^]x - ^^x^ (49) 
cr a(3i) eras (at) 

For p{x, t) > 0, we may set 

p\x,t)^ma^(p\0,t)~-[bit) + 3b{t)^]x-^^x^ o] . (50) 



In conclusion, we have the corresponding functional differential equations (jisp to be the solutions 
of Camassa-Holm equations. 
The proof is completed. ■ 

We notice that the above solutions are not radially symmetric for the density function p with 
b{t) 7^ 0. Thus, the above solutions, cannot be obtained by the separation method of the self-similar 
functional [22], as 

p{x, t) ^ f[^)g{a{U)) and u{x, t) = + b[t). (51) 

On the other hand, for the 2-component Camassa-Holm equations in radial symmetry with 
linear velocity u{r,t): 

Pt + Vpr + pVr = 

(52) 

Vf + WVr + appr = 0, 
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we may replace equation (j22p to have the corresponding step by taking the integration from [0, r] 



^p^ds = -[b{t) + b{t)c{t)] f ds-3[c{t) + c^{t)] f sds. 
2 Jo c's Jq Jq 



(53) 



It is clear for that the rest of proof is similar to have the corresponding result for the solutions in 
radial symmetry: 

Theorem 2 For the 2-component Camassa-Holm equations in radial symmetry there exists a 
family of solutions, 



p\r,t) = max I p\Q,t) - |[6(t) + 36(i)§}]r - -^r\ 



^{r,t)^^J^^r + b{t) 



(54) 



^ A [p2(o,,)] = ^2(0^,) _ + 36(t)^]r - ^^(0,0) = 

where a(3t) and b(t) are the solutions of equations U5]) r and il5\) d. 

3 Blowup or Global Solutions 

To determine that the solutions are global or local only, we can use the corresponding lemma about 
the Emden equation: 



Lemma 3 For the Emden equation I115\) r. 

a(3t) = 

a3 (3t) 

a(0) = ao > 0, d(0) = ai, 
(V ¥ < 0; there exists a finite time T, such that 



(55) 



lim a(M) = 0. 



(2) if ^ — 0, with ai < 0, the solution a{t) blows up in the finite time: 



T = 



ai 



(56) 



(57) 



(3) otherwise, the solution a{t) exists globally. 
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We observe that it is the same lemma for the function a{s) with s = 3i by the separation 
methods in [35]. Therefore, the proofs can be found in Lemma 3 of |22j . 
The gradient of the velocity in solutions (|T5|) and ((54)) . is 

(58) 

When the function a(t) blows up with a finite time T, -g^u{x,T) also blows up at every point x. 
And based on the above lemma about the Emden equation for a(t), it is clear to have the corollary 
below: 

Corollary 4 (la) For ^ < 0, solutions and |5^[ ) blow up in a finite time T; 
(lb) For ^ = 0, with ai < 0, solutions I115\) and blow up in the finite time: 

T=^. (59) 

(2) otherwise, solutions il5\) and |5-^[ ) exist globally. 
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